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Abstract 

We describe in some detail a computer evaluation of a 't Hooft vertex in a two di- 
mensional model using the overlap. The computer result agrees with the known exact 
continuum value, and in this sense our work is a first successful fully dynamical simulation 
of a chiral gauge theory on the lattice. We add some new data to numbers obtained earlier 
and provide a selfcontained description which should make it easy for others to reproduce 
and follow up on our work. 
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1. Introduction 

Unlike for QCD, a communal base of numerical experience about simulating cliiral 
gauge theories on the computer is lacking. Even on what should be done in principle there 
is no consensus; in practice, handling the generically complex action is likely to present 
a formidable problem by itself. One approach for regularizing chiral gauge theories, the 
overlap [1], has been recently subjected to a full dynamical test for a simple model [2]. The 
result was positive and therefore a more detailed description of what was done in practice 
has become relevant. Previous publications emphasized theory and were brief on numerics. 
Here our objective is to collect all the technical details needed to carry out our test, so 
that the reader could reproduce and follow up on our numerical work without having to 
sift through several papers. The reader would be helped by QCD experience only to a 
limited extent; because of the presumed lack of familiarity with the overlap in practice our 
presentation will be quite explicit. 

2. Continuum features. 

Our particular model is a chiral abelian gauge theory in two Euclidean dimensions. 
The theory is defined on a torus of sides h = t, I2 = I- Matter consists of four left-handed 
Weyl fermions, %/, of charge one and one right-handed Weyl fermion of charge two, 
The action is: 

^=4^2/ d'^^^l^ - E / + ^^m)X/ - / d'^xi^a^id^ + 2iA^)i^, (2.1) 

where ui = 1, cr2 = ■i and ji = 1,2. The U{1) gauge symmetry is anomaly free by 
2^ = 1^ -H 1^ -h 1^ + 1^. The boundary conditions are: 

F^i/(^x -\- 1 jj^fi) — F.yi,{x) 

for = 1,2. /i is a unit vector in the jj. direction. The Xf ^^id t/^ fields obey complex 
conjugate boundary conditions. The 6^ are given by: 

6l = 0; 6? = 0; 6?=^; h\=]^; h\ = Q; = ^; 6^ = 0; = ^. (2.3) 

In the t = I = 00 limit the following properties are relevant to the present work: 
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(2.2) 



• The physical Hilbert space can be written as a tensor product of two spaces which are 
not connected by the evolution operator. One of the spaces carries a representation 
of a massless conformal theory and the other describes a massive field theory. If we 
replaced ■^-matter by four right moving Weyl fermions of charge one the theory would 
be vectorial. Again, a similar factorization would hold and the massive factor would 
be isomorphic to the massive factor in our original chiral model. 

• The massless sector consists of six left moving Majorana Weyl fermions forming a 
sextet under the global SU{4:) acting on the /-index of the X/'s- These particles 
are noninteracting. One can choose interpolating fields for these particles which are 
neutral objects local in the original fields: 



P/1/2 
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eo 



(2.4) 



The prefactor is chosen so that p becomes a canonical field at large distances. 7 
is Euler's constant. The mixing on the RHS is possible because fermion number 
is violated in quanta of two by instanton effects. The exact low energy effective 
Lagrangian of the model, written in terms of the p-fields, is 



(2.5) 



/l>/2 



Among the terms on the RHS we have a term proportional to a 't Hooft vertex, V{x)^ 

which we choose to define as: 

.2 



V{x) 



Xi{x)X2{x)xs{x)Xa{x)'4}{x){(7 ■ d)ip{x). 



(2.6) 



This operator has a nonzero expectation value: 

{V)oo = ^ ~ 0.081 



47^3 



On the finite t x I lattice one has instead: 



/T/\ 647r 



{tm^Y 



coth 



'7 



2 ^ 



where = ^ and the functions F{^) and H{^, r) are defined below: 

1 



E 

n>0 



[I- 



Vn2 + (e/27r)2J^ 

H{tT) = y , 



1 



(2.7) 



(2.8) 



(2.9) 
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We shall need {V)txi * = ^ = ^ there its value is 0.0389, substantially smaller 
than the value at infinite volume. 

The objective of our numerical work is to reproduce this last continuum number on 
the lattice. Namely, we wish to evaluate the ratio of formal path integrals: 

inU^dxfdxfmdi;]nU[dA,]eS 

The path integrals are rendered finite by replacing the torus by a toroidal square lattice 
with spacing a. To preserve in variance under rotations by ninety degrees the two sides 
of the lattice are made equal. One wishes to show that the correct continuum number 
emerges in the limit a — > 0. 

The following decomposition of simplifies the path integrals in the continuum: 

27r 27r(7 2tt 

Ai = d2^+ —hi + ig-^dig, A2 = --^xi - dicj) + —h2 + ig-'^d^g. (2.11) 

Above, is a periodic complex valued function on the torus with \g{x)\ = 1, representing 
the gauge degree of freedom. is a real periodic function on the torus and has no zero 
mode: J d'^X(p = 0. hf^ are two real constants restricted to the intervals (—1/2, 1/2] which 
can be thought of as parameterizing two Polyakov loops on the torus, q is an integer 
identifying the topological class. The pure gauge part of S depends only on (f) and q: 



The gauge measure [dA^ is replaced by '^^[dg]d'^h[d(l)\. Because the number of fermionic 
zero modes is determined by g, only the g = 1 term contributes in the sum over q in the 
numerator of (2.10) and only the g = term contributes to the denominator there. Thus, 
the sum over q collapses to a single term. Gauge invariance implies that the integral over 
g can be ignored. Only the fermionic part of the integrand depends on h^. 

3. Lattice formulation 

On the lattice the "measurement" of < F > proceeds by two separate Monte Carlo 
simulations: one for the numerator and the other for the denominator in (2.10). To the 
numerator only gauge fields which induce exactly one fermionic zero mode for each charge 
one fermion and exactly two zero modes for the doubly charged fermion contribute. The 
denominator is entirely determined by configurations with strictly no fermion zero modes. 
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In principle, given a gauge configuration, one needs to perform some calculations to count 
the zero modes. However, by electing to work with a noncompact pure gauge action 
it becomes possible also on the lattice to predict the number of zero modes without a 
calculation simply by setting q. It is still up to the fermionic part of the system to "agree" 
with our prescribed q. We monitor for "disagreements" where the number of fermionic zero 
modes on the lattice is not what it should be for the given q. With our choice of parameters 
such disagreements are so rare that we never encounter them during our simulations. The 
main gain from this is that we can set q by hand and our result does not suffer from the 
statistical noise inherent in simulating the prefactor e~*« from (2.12). 

Our objective is to obtain the expectation value of the 't Hooft vertex in a fixed, 
rather small, physical volume. This is implemented by setting the bare coupling, eo, to 
1.5^ when we work on a finite L x L lattice. The continuum limit is taken by increasing 
L. Note that the parameter was introduced as a bare parameter below equation (2.8). 
A priori, there is no guarantee that the lack of renormalization of the bare coupling in 
the continuum should hold also on the lattice. We assume that it does, and if our final 
result agrees with the continuum value we can say that the numerics are consistent with 
this assumption. 

3 A. Bosonic variables. 
We use the same notation as in the continuum for 0, g, h, q. The pure gauge action is 

where 

A(j){x) = J^{4>{x + jl) + (p{x-jl))-44>{x), 4>{x + Ljl) = (j){x). (3.2) 

All site index components X/^ are between and L — 1 and if they appear outside this 
range in any of our equations one has to use the appropriate boundary conditions to bring 
them back within range. The field (p is constrained by Ylx = 0- Fourier space this 
constraint is trivially implemented. The pure gauge integration measure is taken as: 

Ve"^' / d//--- = Ve"^' fY[dg{x) I ' d^h [ TT #(x)e-^°('^) • • • (3.3) 

and the pure gauge partition function is 
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(3.1) 



Both Monte Carlo simulations carried out to get <V> will calculate the expectation value 
of appropriate "fermion observables" in an ensemble governed by d/i. In other words, both 
the numerator and the denominator of (2.10) are divided by jV. The "fermion observables" 
are functions of g,h,(f) that represent lattice versions of what in continuum is viewed as 
the result of carrying out the Grassmann integral in a fixed gauge background. With q 
fixed the pure gauge action depends only on 0, just like in the continuum. The "fermion 
observable" depends also on g, in addition to /j.^. Hence, the g integral cannot be dropped. 

We generate gauge fields using the pure gauge action and treat the fermion integrals 
in a fixed gauge field background as observables in the pure gauge theory. In our Monte 
Carlo implementation there is a difference between the numerator ond the denominator. 
The measure for generating the variables h is flat in (3.3). We shall keep it this way for 
the numerator. In the denominator however, we replace the flat measure for /i by a mea- 
sure including the known /i-dependent factor from the continuum fermionic contribution. 
This factor is extracted back from the fermionic observable (the equivalent of the chiral 
determinant), so the answer is unchanged. We are thus exploiting the fact that the depen- 
dence on h factorizes from the fermionic determinant in the continuum and from the pure 
gauge action both in the continuum and on the lattice. This trick signiflcantly reduces the 
statistical error of the Monte Carlo simulation "measuring" the denominator. 

For the numerator the h/j, are generated randomly in the interval [—1/2, 1/2]. For the 
denominator the are generated by a weight proportional to: 



1^(2/11, 2/i2)|2 = 



E 



7r(n+2/i2)^+2i7rn2/ii 



(3.5) 



The other gauge variables are simulated in a standard way: Sq defined in (3.1) and 
(3.2) can be written in momentum space as 



1 



L-1 



feo,fei=0 



4sm (— — ) -h 4sm (— — ) 



(3.6) 



where 



4){x) = 



L-1 



E ^(^) 



e ^ 



■k-x. 



(3.7) 



feo,fci=0 

and —k^ = mod(L — fc^, L). If k ^ —k, we generate two random numbers, ui, U2, normally 



distributed according to J duexp 



TV 



4sin2(^) + 4sin2(i^) 



and set (f){k) — ui + 



iu2- If A; = —k, we generate one random number, u, normally distributed according to 
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J duexp 



2tv 



4sin2(^) +4sin2(^)'^ 



several independent (f){x) configurations. 

For the denominator we can therefore replace dfi in (3.3) by 



and set (f){k) — u. We then use (3.7) to obtain 



L/2 J-1/1 ^ 



\d{ri,r2) 



(3.8) 



n#( 



X e 



V <^(.)=oV \0{rur2)\ 



where = 4 /_]^^2 <^^^|^(^i7 ^2)^ and the last factor is thought of as part of the observable. 
In this manner we obtain in the range [—(1/2), (1/2)] distributed according to (3.7). 
Note that q is presumed fixed already. For the numerator, there are no r, Zt variables and 
the factors containing them are absent. Finally, the gauge variables g{x) are U(l) gauge 
transformations randomly generated at each site x. Each gauge configuration we generate 
is statistically independent of the previous ones if we assume a high quality random number 
generator. Wc used RANLUX at luxury level 4. * 

The fermions do not depend on the noncompact fields directly. They only depend 
on ordinary, compact, link variables. Thus, if one wished to replace the noncompact pure 
fermionic action by a compact one the fermionic part of our code would remain unchanged. 
The connection between the variables h, cf), q we generate as described above and the link 
fields transmitted to the fermions is given by: 



Uiix) 



(?(,2-)e^^e['^(^)-'^(^-i)le^^2^*(x + i) if Xi = L - 1 
g{x)e r-^ e^'^^^^~'l'^^~^^^ g* {x + 1) otherwise 

(3.9) 

U2{x) = ^(x)e^e[^(^-2)-'^(^)le-^^^/(a; + 2). 



Note that for g 7^ we take care of the needed twist at the boundary in the first line of 
(3.9). The noncompact character of the pure gauge action is reflected in that the mapping 
{g, h, 0, q} {U^} is many to one. 



* Subtract-and-borrow random number generator proposed by Marsaglia and Zaman, 
implemented by F. James with the name RCARRY in 1991, and later improved by Martin 
Liischer in 1993 to produce "Luxury Pseudorandom Numbers". Fortran 77 coded by F. 
James, 1993 [3]. 
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The right handed fermions hve on the lattice the gauge fields are defined on. The par- 
allel transporters between neighboring sites for these doubly charged fermions are U^{x)'^. 
Up to gauge transformations the U^{x) are quite close to the identity. The U'^ fields will 
be a little farther from the identity. The deviations of z(l — t/^) from 2 A are cutoff effects 
of some typical magnitude. Similar cutoff effects also exist for the singly charged fermions 
and there is no numerical reason to try to make them smaller than the ones for the doubly 
charged fermions. We can thus put the charge one fermions on a twice as coarse lattice, 
thereby saving some computer time without increasing the ultraviolet cutoff effects. 

3B. Fermionic variables. 

The definition of the fermionic observables is in terms of fermionic creation and an- 
nihilation operators associated with the continuum Grassmann variables. Unlike their 
continuum counterparts, these operators are Dirac spinors (not Weyl spinors) thus ending 
up with twice as many components. The charge 2 iIj{x) field is represented by b{x; a) where 
a = 1,2 is the extra spinor index. The charge 1 X/(x) fields are represented by af{y;a). 
The sites y are on the coarse lattice and their integer components range sequentially from 
to — 1. L is always even. Whenever a site index is outside its fundamental range one 
has to use the following boundary conditions to bring it back within range: 

b{x + Lft^-a) ^b{xij,;a), 0/(2/ -/x^; a) = e^''*^a/(2/; a), (3.10) 

where bj^ are given by (2.3). The operators a, b obey canonical anti-commutation relation 
and act on the usual Fock space. and Hy are elementary links on the fine and coarse 
lattices respectively. 

For the numerator we need to discretize (2.6). To this end we introduce the operators 
Vfif2{y) below: 

= «/i(2/;2)a/,(2/;2)6t(2y;2). (3.11) 

The site argument of b is relative to the fine lattice. All operators in (3.11) reside at 
the same physical site. The spinorial indices are chosen to match the continuum Lorentz 
transformation properties of the original fields. Thus, the ?7's are lefthanded and U{1) 
neutral. With the help of the ?7's we can produce a gauge invariant point split discretization 
of V. The r] fields also appear in the definition of the p field in (2.4). 

The central object of the overlap is a Hamiltonian that acts on the Fock space ir- 
reducibly representing the complete algebra of anti-commutators. The Hamiltonian is 
bilinear and separable by flavor. These two properties represent the bilinearity and sep- 
arability of the continuum action. The Hamiltonian (with implicit summations over spin 
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and site indices) is given by: 

4 

7^ = 6tH^(t/'')6-^4H«(t/«)a/. (3.12) 
The single particle Hamiltonians appearing in (3.12) are given below: 

. 2 

C{z, z'- U)^-Y,<y^ - 5.,.'+A. , (3.14) 

-. 4 

B(^, U) = -Y, [2<^..' - <^.',.+A.(t^M(^)) - 5z,z'+[.AU;i^'))] ■ (3.15) 

(3.13) displays the dependence on spin indices in matrix format while (3.14-15) define 
the site dependence explicitly. When using (3.13-15) in (3.12), equation (3.10) must be 
employed to obtain the correct form of 7i, with site index coordinates restricted to [0, L — 1] 
for b and to [0,L/2 — 1] for Uf. U°''^ are the link variables as seen by the appropriate 
fermions. In our case the a fermions reside on the coarse lattice and the U"" link variable 
live on the links of that coarse lattice. The link variables are on the original lattice and 
are chosen as identical to the link variables U introduced in (3.9). The C/" link variables 
are locally defined in terms of the t/'s as follows: The parallel transporter between y and 
y + fiy on the coarse lattice, U^{y)^ is given by: 

W 

W =Ui,{2y)U^{2y + fi^) + U^{2y)U^{2y + Va=)U^{2y + + fia^)U*{2y + 2//^.) 

+ KC^y - ^x)U^{2y - i'j;)Ufj,{2y - + fj,x)Ut,{2y + 2flx - O^), v H- 

(3.16) 

To obtain a proper description for chiral fermions one needs to pick < m < 1. In 
practice, we used m = .5. This is a value empirically tuned to minimize potentially UV 
divergent effects due to a dynamically generated Thirring coupling. 

The fermionic observables are matrix elements of fermionic operators made out of a, b 
in a way corresponding to the continuum expression in terms of Grassmann variables. The 
matrix elements are taken between two different states: The "ket" is the ground state of 
H, |0)j7^, with the phase fixed in a manner which is dependent on the topology: 



WB, 



?(0|(a{(2)ai(l)6(l)6t(2)|0)r (3.17) 
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is real and positive if q = and 



E EE 'M.fsU T(0|(al(2)ai(l)6(l)6t(2)) r^f.fM ^/^(y + A.)|0)r 
/i./2./3./4=i A*— 1 y 

(3.18) 

is real and positive if g = 1. The phase of the state ^^(0| is immaterial since it docs not 
depend on U . The sum over y is over all sites on the coarse lattice. ai{a) — ck) is 

the operator associated with a single particle zero momentum state for the g = 1 fermion 
that has periodic boundary conditions in both directions and the sum on y is over all 
points on the coarse lattice. Similarly, h{cx) — h{x\ a) is the operator associated with 
a single particle zero momentum state for the q — 2 fermion while the sum on x is over 
all points on the fine lattice. These operators create states that carry no momentum. The 
combinations are also charge neutral. All insertions are best thought of as acting to the 
left and defining a U -independent reference state used to fix the phase of the ground state 
of the [/-dependent system. The reference state is a Lorentz scalar and carries no U{1) 
charge. 

The "bra" state in the matrix element representing the result of the Grassmann inte- 
gration, (0|, is [/-independent and given by the ground state of an ultralocal Hamiltonian, 
similar to (3.12) only that the single particle Hamiltonians (3.13) are diagonal having only 
mass terms with the parameter m there replaced by a negative number. The phase choice 
for the "bra" is immaterial as long as it keeps the state [/-independent. 

With these definitions in place we finally can write down the lattice formula for (2.10) 
in the overlap formalism on an L x L lattice as: 

27r2 --2jlL^ N N 



4 2 



^=ki ^ EE^M. f^rnvM-Mmuiy+f^vmr, (3.20) 

/i./2./3./4=i M y 

D=J^ mr- (3.21) 

The sum over y in (3.20) is over all points on the coarse lattice. 

Expressions (3.20) and (3.21) are still not in an explicit enough form to transcribe 
directly into a program. First we need to rewrite (3.20) and (3.21) in first quantized 
forms. To this end, we diagonalize the single particle Hamiltonians (3.13): 

HO = OA. (3.22) 
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The A are diagonal matrices and the real elements along the diagonal are in decreasing 
order top to bottom. We work with five Hamiltonians: four are x matrices and 
differ by the boundary conditions on the unit charged fermions. The fifth is an 2L^ x 2L^ 
matrix representing the doubly charged fermion. 

We also define some gauge field independent reference matrices Oq: These will be used 
in implementing the Wigner-Brillouin phase definition. We shall need a realization of the 
bra state in (3.17) and (3.18) after being acted on by the ai and b operators. 



Ho is one of the five free Hamiltonian obtained by setting the gauge links to unity. The 
boundary conditions are maintained. For a trivial gauge configuration one can diagonal- 
ize the matrices by Fourier transform. Except for fermions which are periodic in both 
directions the ordering of Aq is the same as above. For periodic fermions there are zero 
momentum eigenstates and the action of the di and b operators is implemented by switch- 
ing the position of the positive eigenvalue corresponding to zero total momentum with 
the position of the negative eigenvalue corresponding to zero momentum. This induces a 
switch of two columns in Oq 

It is convenient to partition the O matrices into the following pattern: 



The submatrices have half as many rows as the original O matrices. These rows are 
indexed by a site index, running over values for Oi, (defined on the fine lattice) and 
over values for Ol (defined on the coarse lattice). The first columns of O contain 

the eigenvectors corresponding to positive eigenvalues of H. This set of columns defines 
the horizontal partition of the O matrices. (Exactly zero eigenvalues of H are possible 
but are non-generic.) The Oq reference matrices are partitioned into square submatrices 
OoLL, OoLR, OoRL, OoRR similarly to (3.24). 

We define a configuration to have zero lattice topology (in practice, these are config- 
urations with q = 0) when all submatrices are square. The lattice topological number will 
be unity (in practice, these are configurations with q = 1) when all O^^ are rectangular 
with L^/4 -I- 1 columns and Oll is rectangular with L'^ + 2 columns. When the lattice 
topological number is unity we associate with the rectangular matrices the following site 
dependent square matrices: 



HoOo = OoAo. 



(3.23) 




(3.24) 
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ORR{U^]Xy,H) = {Orr{U^), V{Xy) V{Xy+fly)), 

Olr{U'';0,0) = {Olr{U''), 0, 0). 
The site y is on the coarse lattice. Xy = 2y identifies a site on the fine lattice that also resides 
on the coarse lattice. p,y connects two sites of the coarse lattice. Vj^j^{U°'; y) is the y^^ row 
of 0^j^j^{U°-). v{xy) is a column vector of length that is zero everywhere except at the 
site Xy and "0" denotes a zero column vector. The 0{^{U'';y) are {L^/4: + l) x (L2/4 + 1) 
matrices. ORR{U^\y^ ji) and Olr{U^'-i^,^) are x matrices. The site dependence of 
the square matrices defined in (3.24) corresponds to the sites of the 77 operators in (3.20). 
We need some extra rectangular reference matrices for the case of unit lattice topology: 

<L(0) = (Oif^^, 0), Oi^^{y) ^ {Ol^^, v{y)). (3.26) 

Now we can write down the single particle representations of the matrix elements in 
(3.20) and (3.21). The overlap appearing in (3.21) has the following explicit form: 

mT = \[AeiOi^{U'')d.eiORR{U')e"^^''^ (3.27) 
/ 



with 



In (3.28), 



deiWBRRiU^) 1 T-rp detWB{^(£/") ] 



det WBHfl(t/^) I J y L I det WB{^ (t/«) | -I ' 



V^Brr{U') = [Orr{U')]'Oorr+[Olr{U^)\''Oolr. 



(3.29) 



QRL- 



The overlap appearing in (3.20) becomes (no sum over /j), 

{^\VfuM^hh{y + m)T = ^hhhh^^^O{\{U-;y)d.eiOi\{U'^;y) 
det Oi\{U''- y + fly) det 0{\{U''- y + fiy) det Orr{U''- Xy, i^)e'^^''^ 

with 



(3.30) 



In (3.31), 



^b(«7) = E EE i'hhhff ^MdetWBf^(C/-;y)detWBf^(t/'^;y) 

/l./2./3./4 = l IJ.=1 y 

det WB{\,(C/'^; y + fiy) det WB{''^(t/«; y + Aj/) det WBji,^(t/«; a;^,, //). 



(3.32) 
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The matrices appearing in (3.32) are given by 

WBRR{U'^;Xy,t,) = [ORR{U'^;Xy,fi)]^OoRR+[OLR{U'';0,0)]^OoLR 

3C. The Monte Carlo procedure. 

We have seen that the bosonic variables are integrated over by Monte Carlo while the 
fermionic variables are "integrated" out "exactly" . The bosonic integration has two stages, 
consisting of two Monte Carlo processes, one running within the other. The outer process 
(for either q = or q = 1) averages over gauge orbits, generating the fields (p{x) and h^. 
The links on the coarse and fine lattice arc then defined up to a gauge transformation. 
Setting g{x) = 1 in (3.9) this gives us a representative of an orbit. Loosely, we can refer 
to this as the "Landau" gauge representative, since in the continuum, setting g to unity 
in (2.11) would imply = 0. We then compute the needed determinants (3.27-33), 

depending on the topological sector we arc in. 

The internal Monte Carlo procedure averages over gauge transformations after a "Lan- 
dau" configuration (an orbit) has been selected and the fermionic observables have been 
calculated. For each x we randomly draw a unimodular complex number g{x). When 
the links are replaced by gauge transformed links (see (3.9)) the dependence of the the 
Hamiltonians H"''' in (3.13) on the gauge transformation g can be absorbed into a unitary 
matrix C^'^ which conjugates the "Landau" gauge H"'^. ( ) is a diagonal L^/2 x L^/2 
(2L^ X 2L^) matrix whose rows and columns are indexed by sites and spin components. 
The entries on the diagonal are the appropriate g for the site. Thus the Hamiltonians after 
the gauge transformation are expressed in terms of the "Landau" gauge Hamiltonians as: 

H«'''(gauge transformed links) = [G"'^] H"'''(Landau gauge links) [G"'^. (3.34) 

The effect of the conjugation in (3.34) is to multiply the "Landau" gauge Ol and from 
the left by G"'^. The submatrices (3.24) and (3.25) also get multiplied from the left by 
appropriate diagonal submatrices of G"''^. The determinant factors in (3.27) and (3.30) 
change by gauge factors given by the product of the five determinants of these diagonal 
submatrices. Because of gauge invariance, the dependence on the site y in (3.30) cancels 
out from its gauge factor and both gauge factors turn out to be given by Ylx 9{^)[Yly ^(2/)]^ 
(we remind the reader that the x^s are sites on the fine lattice and the y's are sites on the 
embedded coarse lattice). 

The diagonal matrices implementing the gauge transformations on the matrices O get 
caught in between the matrices multiplied in (3.29) and (3.33) and the related determinants 
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(3.33) 



in (3.28) and (3.32) have to be recomputed. Since no rediagonalization of the Hamiltonians 
is necessary, and since the changes in (3.27) and (3.30) are simple, the total amount of 
computation needed to account for link updates along gauge orbits is smaller than for link 
updates changing gauge orbit. 

Every orbit we generate is completely independent statistically (assuming an ideal 
random number generator) from every other orbit. Similarly, all gauge transformations 
we generate are mutually independent. The absence of correlations implies that statistical 
error analysis is relatively simpler here than in standard Monte Carlo. Of course, since our 
observables are non-local, one needs to look at their distributions to ascertain that they 
are reasonable and that our numerical approach hasn't been rendered hopeless by large 
fluctuations. For sizes L = 8 — 24 the fluctuations are manageable and we can make quite 
accurate statistical determinations. The ability to control the fluctuations is helped by the 
adjustment of the coupling constant with L which keeps the physical volume finite. 

For a given gauge configuration the number of floating point operations is quite high, 
~ lO-^-*^ for the highest L's. The evaluations of the gauge averages of the imaginary parts 
involve some cancelations which are roundoff sensitive. Thus, it is necessary to use double 
precision floating point arithmetic throughout (except in the generation of the random 
numbers) . 



We evaluated the quantities N and D ((3.20-21)) for all even L's between 8 and 24. 
For a single gauge configuration, the evaluation of a D-contribution is faster than the 
evaluation of an A'"-contribution. This is countered somewhat by smaller fluctuations for 
N. We aimed at comparable accuracies in N and D and ended up spending ten times 
as much computer time on evaluating N than on calculating D for the larger L's. For a 
fixed gauge configuration and asymptotic L, N is an L^ computation while D is an L^ 
computation. 

It was determined that reasonable values could be obtained for D with 100 orbits plus 
100 gauge transforms per orbit. For N 20 orbits with 40 gauge transforms sufficed. This 
yielded relative errors (single standard deviation) in N and D of a few percent. 

Numerically, N and D are dominated by ultraviolet effects. When L is large we expect 
the following asymptotics: 



4. Numerical results. 



\og{N) = noL^ + m log L 77,2 + 




(4.1) 



log(L>) = doL^ + di log L + ^2 + 



(is log L (i4 



L2 ' L2 
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no and do are evaluated from the noninter acting theory: 



no = do 



0.242586 . 



(4.2) 



A regularized version cannot reproduce (2.1) exactly in the continuum limit because 
an exactly marginal (in the renormalization group sense) coupling exists and obeys all 
symmetries. This is the "Thirring" coupling: 



gxh can have either sign. Depending on it we can have ni > di or ni < di. 

Generically, by tuning a free parameter in the regularization one can adjust the ap- 
proach to continuum to have rii = di. In our case this is approximatively achieved by 
choosing m = .5 (see paragraph after (3.16)). Of course, the tuning cannot be exact, so 
we will expect N/D to either diverge or go to zero as L becomes truly large. But, as 
long as |ni — (ill logL is numerically negligible log{N/D) is dominated by 712 — d2- The 
first subleading correction is of order 1/L^ since there are no dimension three irrelevant 
operators to contribute to the action or the operator. The exact result quoted in section 
2 is essentially 712 — d2 for the case ni = di. 

Schematically, we would expect the following behavior of the ratio as a function of 




(4.3) 



l/L^: 




Expected Schematic Behavior 



l/L^ 



Figure 1 Effect of Thirring coupling. 
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The data would follow either ABC or ABD, depending on the sign of qth- The better 
we tuned m the closer will B by to E, the exact result for qth = 0. This is only the simplest 
scenario for what happens around B. Since various terms make comparable contributions 
the curves may get more complicated there. 

The corrections that come from the operator (rather than from the action) can 

be evaluated at grh = in the continuum. This yields a correction factor (point-split 
factor) of: 

1*00 7 Q 

ps = exp{-[ —e-*(^+^) + 1.15443133 + 2 log-]}. (4.4) 
Jo ^ L 

The factor ps can be used to improve convergence. The number 1.15443133 above is two 
times Euler's constant and 3/L comes from le^j \pK. We use the complete formula for ps, 
including all subleading terms in Xjl?. 
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Figure 2 Data obtained with gauge averaging. Convergence improves by inclusion 
of the point-split factor ps. 

Figure 2 contains our results for the logarithm of the vertex, log <V >. The exact 
continuum limit is also shown ( log <V > [L = 00) = —3.247 ). Quite clearly, we see 
convergence towards the exact continuum number for L between 8 and 18. ps indeed 
"improves" the observable. The data for L between 20 and 24 show a certain amount of 
instability. It is conceivable that this reflects region B in our schematic figure. We cannot 
say anything conclusive without better statistics and additional results for higher L's. 
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Figure 3 Data obtained using the absolute value and no gauge averaging. For L = 
18, 20, 22, 24 higher statistics points are added. 

Most of the computation time is spent on gauge averaging. By discarding gauge 
averaging we can increase statistics. A simple gauge invariant way to avoid gauge averaging 
is to replace the integrands in (3.20) and (3.21) by their absolute values. The corresponding 
quantities in the continuum are positive. In figure 3 we show data obtained using the 
absolute values exploiting the same runs as in figure 1. We also add four more data points 
in the range L = 18 to 24 with higher statistics: is evaluated from 100 orbits and D from 
200. We see that the higher statistics, while supportive of some structure around L = 20, 
still do not lead us to a definite confirmation of the qualitative scenario that would relate 
it to a Thirring term. 

To be sure that using the absolute value rather than the gauge averaged results is not a 
source of contamination as far as the possibility of a structure around L = 20 is concerned 
we compare in figure 4 the two (absolute value and gauge averaged) vertices for the entire 
set of available orbits. We see that the difference is quite insignificant numerically and, 
therefore, how gauge invariance is recovered is unimportant. 

In figure 5 we display N and D individually. We subtract the large numbers uqL'^ = 
cIqL^ from the logarithms of N and D and plot the results together with the logarithm of 
the ratio {log{N/D) = 0.474 in the continuum limit). We see that the bulk of D and N is 
given by the free field terms we subtracted, and that the remainder consists of numbers of 
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Figure 4 A comparison of data obtained with gauge averaging to data obtained with- 
out for the same sample of orbits. 
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Figure 5 log N and log D after subtracting the quadratic divergence are small num- 
bers. The logarithm of N/D is also shown. The numbers were obtained 
with gauge averaging. 
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order unity. 

5. Summary. 

Although the tuning of the Thirring couphng is not fuUy understood numericaUy, 
overaU, our simulation confirms the validity of the overlap regularization method of the 
model in (2.1). More work would be needed to understand the possible role of a dynamically 
generated qth- Other two dimensional models can be tested by similar methods. Going 
to higher dimensions would require some algorithmical improvements, as the simulations 
presented here took a considerable amount of time on machines of various types. 

Based on our work, we express our hope that the problem of simulating chiral gauge 
theories has been turned from one of principle to one of practicality. 
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